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Abstract 

The growth rate of matter perturbation and the expansion rate of the Universe can be used to 
distinguish modified gravity and dark energy models in explaining the cosmic acceleration. The 
growth rate is parametrized by the growth index 7. We discuss the dependence of 7 on the matter 
energy density f2 and its current value Oo for a more accurate approximation of the growth factor. 
The observational data, including the data of the growth rate, are used to fit different models. 
The data strongly disfavor the Dvali-Gabadadze-Porrati model. For the dark energy model with 
a constant equation of state, we find that Qq = 0.27 db 0.02 and w = —0.97 db 0.09. For the 
ACDM model, we find that 7 = 0.64±^. For the Dvali-Gabadadze-Porrati model, we find that 
7 — u.oo_ 013 . 
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I. INTRODUCTION 



The discovery of late time cosmic acceleration [l| challenges our understanding of the 
standard models of gravity and particle physics. Within the framework of Friedmann- 
Robertson- Walker cosmology, an exotic energy ingredient with negative pressure, dubbed 
dark energy, is invoked to explain the observed accelerated expansion of the Universe. One 
simple candidate of dark energy which is consistent with current observations is the cos- 
mo logical constant. Because of the many orders of magnitude discrepancy between the 
theoretical predication and the observation of vacuum energy, other dynamical dark en- 
ergy models were proposed [2j . By choosing a suitable equation of state w = p/ p for dark 
energy, we can recover the observed expansion rate H(z) and the luminosity distance red- 
shift relation di{z). Current observations are unable to distinguish many different dark 
energy models which give the same (Il(z), and the nature of dark energy is still a mystery. 



Many parametric and nonparametric mode 
the property of dark energy 



-independent methods were proposed to study 
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Bear in mind that the only observable effect of dark energy is through gravitational 
interaction; it is also possible that the accelerated expansion is caused by modification of 
gravitation. One example of the alternative approach is provided by the Dvali-Gabadadze- 
Porrati (DGP) brane-world model [3], in which gravity appears four dimensional at short 
distances but is modified at large distances. The question one faces is how to distinguish such 
an alternative approach from the one involving dark energy. One may answer the question 
by seeking a more accurate observation of the cosmic expansion history, but this will not 
break the degeneracies between different approaches of explaining the cosmic acceleration. 
Recently, it was proposed to use the growth rate of large scale in the Universe to distinguish 
the effect of modified gravity from that of dark energy. While different models give the same 



late time accelerated expansion, the growth of matter perturbation they produce differ [19]. 
To linear order of perturbation, at large scales the matter density perturbation 5 = 5p/p 
satisfies the following equation: 

5 + 2H5 -4irG effP 5 = 0, (1) 



where p is the matter energy density and G e ff denotes the effect of modified gravity. For 
example, G e ff/G = (4 + 2u;)/(3 + 2uS) for the Brans-Dicke theory [20j and G e ff/G = 
2(1 + 2fi 2 )/3(l + VI 2 ) for the DGP model Q, the dimensionless matter energy density 



Q = 87cGp/(3H 2 ). In terms of the growth factor / = din 5 /din a, the matter density 
perturbation Eq. ([[]) becomes 



n. 



(2) 



where /' = df/dlna. In general, there is no analytical solution to Eq. (121) , and we need 
to solve Eq. (TS]) numerically; it is very interesting that the solution of the equation can 



be approximated as / = Q"' 
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251 ] and the growth index 7 can be obtained for 



some general models. The approximation was first proposed by Peebles for the matter 
dominated universe as f(z = 0) = oIq 6 [22I; then a more accurate approximation, f(z = 
0) = 1 for the same model was derived in 23, 24]. For the Universe with a cosmo logical 
constant, the approximation f(z = 0) = Qq 6 + fi\o(l + ^o/2)/70 can be made [26(. For a 
dynamical dark energy model with slowly varying w and zero curvature, the approximation 
f(z) = f2(2i) 7 was given in 25J. For the DGP model, 7 = 11/16 27|. Therefore, instead 
of looking for the growth factor by numerically solving Eq. (j2J), the growth index 7 may 
be used as the signature of modified gravity and dark energy models. It was found that 
7 = 0.55 + 0.05[1 + w(z = 1)] with w > -1 and 7 = 0.55 + 0.02[1 + w(z = 1)] with 
w < — 1 for the dynamical dark energy model in flat space 27], [28]. Recently, the use of the 
growth rate of matter perturbation in addition to the expansion history of the Universe to 



differentiate dark energy models and modified gravity attracted much attention 
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The dependence of 7 on the equation of state w has received much attention in the 
literature; we discuss the dependence of 7 on Q and Qq for a more accurate approximation 
in this paper. We discuss a more accurate approximation for the dark energy model with 
constant w in Sec. II. Then we discuss the DGP model in Sec. III. In Sec. IV, we apply 
the Union compilation of type la supernovae (SNe) data Q, the barvon acoustic oscillation 
(BAO) measurement from the Sloan Digital Sky Survey (SDSS) 46], the shift parameter 



measured from the Wilkinson Microwave Anisotropy Probe 5 yr data (WMAP5) 



Hubble parameter data H(z) [48|, |49j, and the growth factor data f(z) [42l . 



44 



ijtthe 
5o| to 



constrain the models. We also use the growth factor data alone to find out the constraint 
on the growth index 7. We conclude the paper in Sec. V. 
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II. ACDM MODEL 



We first review the derivation of 7 given in 25j. For the flat dark energy model with 
constant equation of state w, we have 



— |[i+«d 

The energy conservation equation tells us that 



n)]. 



ft' = 3wft(l - ft). 



Substituting Eqs. flUJ and gj) into Eq. (J2D, 

3 w ft(l-ft)|- + / 2 + 
ail 

Plugging / = ft 7 into Eq. (jSJ), we get 

g?7 



we get 
1 3 



w(l-ft) 



f 



-ft. 



3wft(l -ft) In ft 



dft 



3wft( 7 - 1/2) + ft 7 - ^ft 1 " 7 + 3w T - + X - 



Expanding Eq. around ft = 1, to the first order of (1 — ft), we get 25] 

3(1 -w) 3 (1 -w)(l - 3w/2) 



7 



(3) 



(4) 



(5) 



(6) 



(7) 



5-Qw 125 (1 - Qw/5) 2 
In order to see how well the approximation ft 7 fits the growth factor /, we need to solve Eq. 
([6]) numerically with the expression of ft. The dimensionless matter density is 

fto 



ft 



(8) 



ft + (l-fto)(l + ^) 3w 

Since 7 does not change very much, we first use 7 = 700 = 3(1 — w)J (5 — 6w) to approximate 
the growth factor. For convenience, we choose fto = 0.27, and the result is shown in Fig. [U 
From Fig. (TJ we see that ft 7 °° approximates / better than 2%. 
For the ACDM model, w — — 1, so the growth index becomes 



3(1 -w) 



7o 



6 

11' 



and 



7i 



g?7 
rift 



5 — 6w 
3 (l-w)(l-3w/2) 



125 (1 - 

For the ACDM model, the matter density is 

fto(l 



15 
1331' 



ft 



ft (l + 2) 3 + l-ft c 



(9) 



(10) 



(11) 




FIG. 1: The relative difference between the growth factor / and 7o ° with 7^ = 3(1 — w)(5 — Qw) 
for the dark energy model with constant w. 

Substituting Eq. (fTTj) into Eq. (jfJJ) and solving the equation numerically, we compare the 
numerical result / with the analytical approximation fi 7 . The results are shown in Figs. [2] 
and [31 Since 71 is very small, 7 does not change much, and we can take the approximation 
7 = loo- In Fig. [2], we compare f2 7 °° with / for different values of f2 . From Fig. [21 we see 
that f2 7o ° approximates / very well; the smaller Q , the larger the error. When the Universe 
deviates farther from the matter dominated era, the error becomes larger. For f2 = 0.2, 
the approximation overestimates / by only 2%, or 700 underestimates 7. To get a better 
approximation, we expand to the first order of (1 — Q) and use 7 = 7^ + 71(1 — ft). In 
Fig. [21 we plot the relative difference between f2 7 and /. From Fig. [31 we see that using 
7 = 7oo +7i(l — ^) approximates the growth factor much better; now the error is only 0.6% 
for = 0.2. 
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FIG. 2: The relative difference between the growth factor / and W°° with 7^ = 6/11. 
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FIG. 3: The relative difference between the growth factor / and f2 7 with 7 = 700 +71(1 — Q) for 
the ACDM model. 
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III. DGP MODEL 



For the flat DGP model, we have 

G eff _ 2(1 + 20 2 



G 3(1 + 2 )' 

The Friedmann equation tells us that 

H 30 

1p ~ ~TTo' 

The energy conservation equation tells us that 



The matter energy density is given by 



Substituting Eqs. (pL (jlgp and (jMj) into Eq. (J2J), we get 



1 + dO J 1 + 0' l + O 2 
Plugging / = 7 into Eq. (fill , we get 



(12) 



(13) 



t y = - 3n ( 1 -"). (i4) 

1 + v ; 



n= ^o(l + ^) 3 (15) 

[(1 - O )/2 + v/O (l + + (1 - O ) 2 /4] 2 ' 



30(1 -Q)df , 2 2-0 0(1 + 20 2 ) 



30(1-0) In O rf 7 _ 3(1-0)7 2-0 _ 0^(1 + 20 2 ) = 

1 + dO 1 + 1 + l + O 2 ' lj 

Expanding Eq. (1171) around O = 1, to the first order of (1 — 0), we get 

7 =— + —(1-0). (18) 
' 16 5632 v ; v ; 

So 7oo = 11/16 and 71 = 7/5632. The change of 7 is very small because 71 is very small; 
we first approximate 7 by 7oo and the result is shown in Fig. HI From Fig HI we see that 
the error becomes larger when the Universe deviates farther from the matter domination. 
When O = 0.2, 7o ° underestimates / by 4.6%. So the growth factor is overestimated. If 
we use the first order approximation, the error becomes larger because 71 > 0. Linder and 
Cahn give the approximation 27j] 

7 + 50 + 70 2 + 30 3 11 7 M ^ 
7 " (1 + 0^(11 + 50) "l6 + 256 (1 -^ (19) 
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Again, to the first order approximation of (1 — O), the error becomes larger than that with 



7c 



In 



431 ]. the author found the approximation 

11 3 
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;i-n) 
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(20) 



16 lQ y ' 16 

In this approximation, the correction to is too big because 71 is too large even though 
the sign is correct. To get better than 1% approximation, we first assume 7 is a constant 
and solve Eq. (fTTj) at Q = Q to get 7 ; we then approximate 71 = (70 — 7oo)/2(l — fl ). 
The values of 70 and 71 for different Qq are listed in table [I] The difference between with 
7 = loo + 7i(l — ^) and / is shown in Fig. As promised, the error is under 1%. 
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FIG. 4: The relative difference between the growth factor / and J7 7o ° with j c 
DGP model. 
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IV. OBSERVATIONAL CONSTRAINTS 

Now we use the observational data to fit the dark energy model with constant w and the 
DGP model. The parameters in the models are determined by minimizing \ 2 — Xsn + X" 'f 



Xr + Xh + Xf- For the SNe data, we use the reduced Union compilation of 307 SNe [45 ]. 
The SNe compilation includes the Supernova Legacy Survey [51j and the ESSENCE Survey 
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FIG. 5: The relative difference between the growth factor / and IP with 7 = +71(1 — £1) for 
the DGP model. 
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7o 


7i 


0.2 


0.639 


-0.0305 


0.27 


0.648 


-0.0272 


0.3 


0.652 


-0.0257 


0.4 


0.663 


-0.0205 



TABLE I: The values of 70 and 71 in the DGP model for different Oq- 



521 . |53| |. the older observed SNe data, and the extended dataset of distant SNe observed 
with the Hubble space telescope. To fit the SNe data, we define 



X 2 

Asn 



307 

£ 

i=l 



(21) 



where the extinction-corrected distance modulus fi(z) = 5 log 10 [(iL(^)/Mpc] + 25, a. t is the 
total uncertainty in the SNe data, and the luminosity distance is 

f z dz' 

d L (z) = H \l + z) j( (22) 

The dimensionless Hubble parameter E(z) = H(z)/H = [Q {l+z) 3 + (l-Q )(l+z)^ 1+w ^ 2 
for the dark energy model with constant w and E(z) = [f2 (l+z) 3 + (l— fi ) 2 /4] 1 ^ 2 + (l— fi )/2 



for the DGP model. The nuisance parameter Hq is marginalized over using a flat prior. 



To use the BAO measurement from the SDSS data, we define x\ 
0.469(0.96/0.98)- a35 ) 2 /0.0 1 7 2 , where the distance parameter Q, Q 

1/3 

= 0.469(0.96/0.98)-°' 35 ± 0.017. 



A 



0.35 



0.35 



£(0.35) 



0.35 



dz 



E z) 



To use the shift parameter measured 
1.715) 2 /0.021 2 , where the shift parameter 



rom the WMAP5 data, we define x 



3 



n 



1089 



dz 



1.715 ±0.021. 



(A - 

(23) 
(K- 

(24) 



Simon, Verde, and Jimenez obtained the Hubble parameter H(z) at nine different red 



481 ] . Recently, the authors in 



shifts from the differential ages of passively evolving galaxies 

4<J obtained H(z = 0.24) = 83.2 ±2.1 and H(z = 0.43) = 90.3 ± 2.5 by taking the BAO 
scale as a standard ruler in the radial direction. To use these 11 H(z) data, we define 



Xh 



11 

E 



[#<*.(*) - H(z,)}- 



i=l 



a 



hi 



where a hi is the la uncertainty in t 
km/s/Mpc given by Freedman et al. 



re H(z) data. We also add the prior Hq 



(25) 
72 ± 8 



54]. 



For the growth factor data, we define 
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E 

i=l 



[f obs { Zi ) - w{zi)Y 



a 



(26) 



fi 



where 07 ,; is the la uncertainty in the f(z) data. For reference, we compile the available 



data 



42 



n 



|44j, [50|] in Table [TTl The data are obtained from the measurement of the redshift 
distortion parameter (3 = f/b, where the bias factor b measures how closely galaxies trace 
the mass density field. Note that some of the measured redshift distortion parameter (3 
is obtained by fitting f2 06 with Q given by the ACDM model, and some analyses tried to 



account for extra distortions due to the geometric Alcock-Paczynski effect [501 ] . With these 
caveats in mind, it is still worthwhile to apply the data to fit the models. As discussed in the 
previous sections, we can use 7 = 700 within the accuracy of a few percent. For the ACDM 
model, we use 7 = 6/11 and Q given by Eq. ( JTT1) . For the DGP model, we use 7 = 11/16 
and fl given by Eq. (|T5|) . 
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z 


fobs 


References 


0.15 


0.49 ±0.1 


[50, 551 


0.35 


0.7 ±0.18 


[56] 


0.55 


0.75 ±0.18 


[57] 


0.77 


0.91 ±0.36 


[50] 


1.4 


0.9 ±0.24 


[58] 


3.0 


1.46 ±0.29 


[59] 


2.125 - 2.72 


0.74 ±0.24 


[60] 


2.2 - 3 


0.99 ± 1.16 


[61] 


2.4-3.2 


1.13 ± 1.07 


[61] 


2.6-3.4 


1.66 ± 1.35 


[61] 


2.8-3.6 


1.43 ± 1.34 


[61] 


3-3.8 


1.3 ± 1.5 


Liiil 



TABLE II: The summary of the observational data on the growth factor /. 

By fitting the dark energy model with constant w to the combined data, we get \ 2 = 
325.17, Q = 0.2721^1, and w = -0.97 ± 0.09. The la, 2a, and 3a contours of tt and 
w are shown in Fig. El From Fig. El we see that the ACDM model is consistent with the 
current observation. By fitting the ACDM model to the combined data, we get x 2 = 325.48 
and tt = 0.273±0.015. By fitting the DGP model to the combined data, we get \ 2 = 350.12 
and tt = 0.278 ± 0.015. 

If we fit f2 7 to the growth factor data f(z) alone, we can get a constraint on the growth 
index 7. For the ACDM model with the best fit value flo — 0.273, we find that x 2 — 4.55 and 
7a = 0.641q}5. The theoretical value 700 = 6/11 = 0.55 is consistent with the observation at 
the la level. For the DGP model with the best fit value flo — 0.278, we find \ 2 = 5.68 and 
Idgp — 0.55lgj3 which is barely consistent with the theoretical value 700 = 11/16 = 0.6875 
at the la level. The best fit curves W and the observational data are shown in Fig. [71 
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FIG. 6: The la, 2a, and 3a contours of f2 and w by fitting the dark energy model with constant 
w to the combined data. The point with + denotes the best fit value. 

V. DISCUSSIONS 

The simple analytical formulas f2 7 can be used to approximate the growth rate /. The 
value of 7 provides useful information about the dark energy model and the modification of 
gravity. For the dark energy model within general relativity, 7 ~ 0.55. For the DGP model, 
7 = 0.6875. If the accuracy of the growth factor data is in the range of a few percent, we 
can use a constant 700 to approximate /. For example, 700 = 3(1 —w)/(5 — 6w) for constant 



w 



25j . 7oo = 0.55 + 0.05[1 + w(z = 1)] for dynamical models [28J], and 700 = 11/16 for 
the DGP model. The value of 700 is obtained by approximating the differential equation 
of / around Q — 1. This approximation is reasonably good at high redshift (z ^ 1) when 
n « 1. However, at lower redshift the dark energy or the effect of modified gravity starts to 
dominate and Q deviates from 1; we expect the approximation to break down. Therefore, 
to get a better than 1% fit, the Q dependence of 7 needs to be considered. We show 
that 7 = 7oo + 7i(l — fi) can approximate / to better than 1%. The value of 71 is very 
small compared with 7^ and usually depends on the value of fl . For the ACDM model, 
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FIG. 7: The growth factor data and the best fit curves. The solid line is for the dark energy model 
with n = 0.272 and w = -0.97. The dashed line is for the ACDM model with f} = 0.273 and 
7 = 0.64. The dotted line is for the DGP model with £l = 0.278 and 7 = 0.55. 

71 = 15/1331. For the DGP model, we give a prescription of how to find 71; the values of 
71 are listed in Tabled for some values of fi . 

To distinguish different dark energy models and modified gravity, we use the observational 
data to fit the models. Fitting the combined SNe, SDSS, WMAP5, H(z), and f(z) data to 
the dark energy model with constant w, we find that x 2 = 325.17, fio = 0.272io;o22? an d 
w = -0.97 ±0.09. For the ACDM model, we find that y 2 = 325.48 and Q = 0.273 ±0.015. 
For the DGP model, we find that y 2 = 350.12 and fi = 0.278 ± 0.015. The results suggest 
that the data strongly disfavor the DGP model. If we fit the SNe data alone to the DGP 
model, we get y 2 = 313.27 and f2 = 0.186. This result is consistent with the results in 



62 



63 



64 



65 



66l . 1671 . |68j. If the same data were fitted to the ACDM model, we would 
get x 2 = 312.18. If we fit the ACDM model to the combined SNe and H(z) data, we get 
X 2 = 319.45. For the DGP model, we get y 2 = 321.31 and f2 = 0.194. These results show 
that both the ACDM model and the DGP model give almost the same expansion rate. If 
we fit the ACDM model to the combined SNe, H(z), and f(z) data, we get x 2 = 324.53. 
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For the DGP model, we get \ 2 = 331.53 and fi = 0.207. With the addition of f(z) data to 
the expansion data, the DGP model is readily distinguishable from the ACDM model. As 
mentioned above, the best fit value of Qq in the DGP model tends to be lower, i.e., Qq ~ 0.2, 
when we fit the model to the SNe, H(z), and f(z) data. On the other hand, we get x 2 — 1-4 
and Qq = 0.33 if we fit the model to the distance parameter A and the shift parameter 1Z. 
If the two sets of data are combined together, Qq takes the value in the middle and we get 
a large value of \ 2 ■ When we fit the models to the combined SNe, the distance parameter 
A, the shift parameter 1Z, and H(z) data, we get Q = 0.274 and \ 2 — 320.6 for the ACDM 
model and Qq — 0.275 and y 2 = 343.5 for the DGP model. This result is consistent with 
the analysis by Song et al 62j. In 62|, they found that the flat DGP model is excluded at 
about 3a by the combined SNe, 3 yr WMAP, and Hubble constant data. 

The observational data of f(z) can be used to provide information on the growth index 7 
and the modified gravity. As discussed above, 7 is almost a constant; we can use W with a 
constant 7 to approximate f(z). For the ACDM model, we find that 7a = 0.64^^5, which is 



consistent with the theoretical value 0.55. This result is also consistent with that in |42l. 144 ] . 
For the DGP model, we find that 'Jdgp — 0-55l ' 13 . The theoretical value 7 = 0.6875 lies 
on the upper limit of the la error. From Fig. [TJ we see that the growth rates for the DGP 
model and the ACDM model are distinguishable even with the best fitting values of 7. If 
the theoretical values of 7 are used, the difference will be larger. The approximation of W 
to / is good for theories with / < 1 and Q < 1. If the observational data for f(z) are 
larger than 1 at high redshift (z ^ 1), then the approximation f2 7 is broken and the effect 
of modified gravity is explicit. For the Brans-Dicke theory, during the matter domination, 
VI = (3lu + 4)(2lo + 3)/6(l + uof > 1 and / = (2 + w)/(l + u) > 1 [|2|; here u is the 
Brans-Dicke constant. In this case, modification of the approximation is needed. This will 
be discussed in future work. In conclusion, more precise future data on f(z) along with the 
SNe data will differentiate dark energy models from modified gravity. 
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